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ABSTRACT 

Let  V  be  (the  boundary  of)  a  simple  polygon  with  n  vertices.  For  a 
vertex  p  oiV,  let  <f>{p)  be  the  set  of  points  on  V  that  are  farthest  from  p, 
where  the  distance  between  two  points  is  the  length  of  the  (Euclidean) 
shortest  path  that  connects  them  without  intersecting  the  interior  of  V. 
In  this  paper,  we  present  an  0(n  log  n)  algorithm  to  compute  a  member 
of  <f>{p)  for  every  vertex  p  of  V.  As  a  corollary,  the  external  diameter  of 
V  can  also  be  computed  in  the  same  time. 
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Introduction 


1      Introduction 

Shortest  path  planning  in  an  obstacle-cluttered  environment  is  a  classical  problem  in 
computational  geometry.  In  two  dimensions,  it  is  common  to  model  the  obstacles  by 
closed  simple  polygons  Vi,  . . .  ,  Vm  whose  interiors  represent  forbidden  regions,  and 
to  model  a  robot  by  a  point.  A  feasible  path  of  the  robot  in  such  an  environment 
is  a  path  that  does  not  meet  the  interior  of  the  obstacles.  Given  two  points  in 
the  plane,  the  distance  between  them  is  the  length  of  the  shortest  feasible  path 
connecting  them;  if  either  of  the  points  lies  inside  an  obstacle,  this  distance  is 
assumed  to  be  infinite. 

The  set  Ui  ^t  is  called  the  obstacle  space.  A  useful  measure  of  the  obstacle 
space  is  the  maximum  distance  between  any  two  points  lying  on  its  boundary.  This 
maximum  distance  is  called  the  diameter  of  the  obstacle  space.  In  this  paper,  we 
axe  concerned  with  the  case  of  a  single  obstacle. 

Given  a  simple  n-gon  V  and  two  points  p  and  q  on  its  boundary,  we  define  the 
external  distance  between  p  and  q  to  be  the  length  of  a  shortest  path  that  joins 
p  and  q  without  intersecting  the  interior  of  V.  For  a  vertex  p  of  V,  let  (f)(p)  be 
the  set  of  points  on  (the  boundary  of)  V  that  are  farthest  from  p  with  respect  to 
the  external  distance,  and  let  the  external  diameter  of  V  be  the  maximum  external 
distance  between  any  two  points  of  V.  The  problem  of  computing  the  external 
diameter  was  first  considered  by  Samuel  and  Toussaint  [ST87].  They  presented  an 
0{n^)  time  algorithm  to  solve  the  problem,  where  n  denotes  the  number  of  vertices 
of  P. 

In  this  paper  we  provide  an  0{n  log  n)  time  and  0{n)  space  algorithm  for  com- 
puting a  point  of  <f){p)  for  every  vertex  p  oi  V.  A  result  of  Samuel  and  Toussaint 
[ST87]  then  shows  that  the  external  diameter  of  V  can  also  be  computed  in  time 
0(n  log  n). 

This  paper  is  divided  into  five  sections.  Section  2  discusses  the  basic  geometric 
concepts  that  we  use  in  this  paper.  In  Section  3  we  prove  some  properties  of 
external  shortest  paths,  which  lead  us  to  an  efficient  algorithm  for  computing  the 
external  farthest  neighbors  for  every  vertex  of  the  polygon.  Section  4  describes 
an  0{n\ogn)  algorithm  to  compute  external  farthest  neighbors.  We  conclude  with 
some  final  remarks  in  Section  5. 


2      Geometric  preliminaries 

Let  V  be  (the  boundary  of)  a  simple  polygon  with  n  vertices.  Let  po,pi,  . . .  ,Pn-i 
denote  the  vertices^  of  V  ordered  in  clockwise  direction  around  it.  For  distinct 
points  p,q  E  V,  let  Vlp,q]  be  the  section  of  V  clockwise  between  p  and  q  inclusive. 


^Indices  of  p,  are  modulo  n. 
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7I"l(p,9)  '\  /'  7rfl(p,9) 


Figure  1:  Left  and  right  shortest  paths 

Define  V{p,q)  =  V\p,q]  -  {p,g};  define  V\p,q)  and  'P{p,q]  similarly.  Let  CHiV) 
denote  the  boundary  of  convex  hull  of  V,  and  qo,q\,  ...  ,qm-i  the  vertices^  of  V, 
lying  on  CHiV)  in  clockwise  order.  Without  loss  of  generality  we  can  assume  that 
qo  =  Po-  'P\qi-,qi+i)  is  denoted  by  K,  and,  \iV[q,,q,j^\\  ^  q^q^+x,  the  closed  polygonal 
region  bounded  by  K-i  and  5,5,+!  is  denoted  by  AC,.  If  q,q,+i  =  V[qi,qi+i],  then  let 
AC,  =  q,qi+i.  Somewhat  abusing  the  notation,  we  will  refer  to  either  AC,  or  AC,  as 
a  pocket.  Note  that  pockets  ACq,  fCi,  ...  ,Km-i  form  a  partition  of  V.  An  edge 
5,5,+i  C  C1-L{V)  is  called  the  lid  of  a  pocket  and  the  points  g,  and  g,+i,  lid  vertices 
of  AC,  (or  of  ACi). 

For  a  pair  of  points  p,  q  not  lying  in  the  interior  of  V,  a  path  from  p  to  g  is  an 
external  path  if  it  avoids  the  interior  of  V.  The  external  distance  d{p,  q)  between 
two  points  p  and  q  is  the  length  of  a  shortest  external  path  connecting  p  to  9  (that 
such  a  path  always  exists  can  be  shown  by  an  easy  compactness  argument).  Since 
the  exterior  of  V  in  the  plane  is  not  simply  connected,  there  are  infinitely  many 
homotopy  classes  of  paths  connecting  two  points  p  and  q  on  V.  Informally,  the 
homotopy  class  of  such  a  path  determines  how  many  times  it  winds  around  V  and 
the  direction  in  which  it  winds  around  V.  It  is  not  difficult  to  show  that  any  path 
connecting  p  e  V  to  q  e  V  that  winds  around  V  more  than  one  complete  revolution 
must  self-intersect  and  thus  cannot  be  a  (globally)  shortest  external  path  between 
p  and  q  ^  p.  We  will  therefore  restrict  our  attention  to  paths  homotopic  to  Vlp,q] 
and  V[q,p],  referring  to  them  as  left  and  right  paths  from  p  to  5,  respectively.  If  p 
=  q,  then  notion  of  left  and  right  paths  is  undefined.  For  an  external  path  u,  let 
|a;|  denote  its  length  in  L2  metric.  For  distinct  p,q  £  V,  we  define  niip^q)  (resp. 
7rR(p,9))  to  be  the  shortest  left  (resp.  right)  path  from  p  to  q.  (Lemma  3.2  below 
shows  that  there  is  a  unique  shortest  left  (resp.  right)  path  from  p  to  q.)  Note  that 
7!'l(p,9)  =  7rfl(g,p),  for  any  pair  of  distinct  points  p,q  eV.  Let  di{p,q)  =  |7ri(p,g)| 
and  dR{p,q)  =  \T^R{p,q)\.    As  previous  discussion  indicates,  there  are  at  most  two 


^Indices  of  qj  are  modulo  m. 

External  Diameter  April  I    1989 


Geometric  preliminaries 


distinct  external  shortest  paths  between  two  distinct  points  p,  g  €  "P,  namely  ^^{p,  q) 
and  T^nipiq)  (see  Figure  1).  In  particular, 

d{p,q)  =  m\n{dL{p,q),dR{p,q)]. 

U  d lip,  q)  <  dR{p,q),  define  Tr{p,q)  =  ttlCp,?),  otherwise  7r(p,g)  =  '^r(p^q)-  Define 
n{p,p)  as  the  "null  path"  from  p  to  p. 

For  a  given  point  p  G  "P,  an  external  farthest  neighbor  of  p  is  a  point  q  E  V  such 
that 

d(p,q)  =  sup{d{p,p')\p  e  V}. 

Let  (f>{p)  denote  the  set  of  all  external  farthest  neighbors  of  p.  We  define  a  more 
specialized  notion  of  external  farthest  neighbors.  Let  (f>'{p)  denote  the  set  of  the 
farthest  neighbors  of  p  among  points  not  lying  in  the  same  pocket  as  p.  Formally, 
if  p  G  /C,,  then 

<t>'ip)  =  {x  G  P[g,+i,9>]  |Vy  G  V[qi+i,qi],  d{p,x)  >  d{p,y)}. 

The  external  diameter  of  V,  denoted  by  T>{V)^  is 

V(V)  =  snp{d{p,q)\p,qeV}. 

A  pair  of  points  p,q  £  V  with  d{p, q)  —  T>{V)  is  referred  to  as  a  diametric  pair.  A 
result  of  Samuel  and  Toussaint  [STS7]  shows  that 

Lemma  2.1  ([ST87])  For  any  given  simple  polygon  V,  there  exists  a  diametric 
pair  p,q  €.  V  such  that  p  (but  not  necessarily  q)  is  a  vertex  of  V. 

D 

Remarks: 

(i)  It  is  easy  to  see  that  there  exist  simple  polygons  in  which  the  external  farthest 
point  of  every  vertex  is  an  interior  point  of  an  edge,  e.g.  consider  any  regular 
k-gon  with  odd  number  of  vertices  (see  Figure  1). 

(ii)  There  are  simple  polygons  in  which  all  diametric  pair  of  points  lie  within  a 
pocket  (see  Figure  2) 

Lemma  2.1  implies  that  a  fast  algorithm  to  compute  external  farthest  neighbors 
immediately  yields  an  efficient  procedure  for  computing  the  external  diameter. 
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Figure  2:  External  diameter  is  realized  by  a  pair  of  points  p,  q  lying  in  the  same 
pocket 

3      Properties  of  external  shortest  paths 

In  this  section  we  study  external  shortest  paths  and  prove  some  of  their  properties 
which  lead  us  to  an  efficient  algorithm  for  finding  the  farthest  neighbors.  First  we 
state  without  proof  some  simple  observations: 

•  The  external  distance  d(p,r)  is  a  metric.  In  particular,  d{p,r)  +  d{r,s)  > 
d{p,  s)  for  any  p,r,s  G  V. 

•  For  any  two  points  p,r  £  V,  n{p,r)  is  a  polygonal  chain  which  lies  entirely 
within  the  region  enclosed  by  C7i{V),  and  whose  vertices  occur  only  at  vertices 
oiV. 

•  An  external  shortest  path  between  two  points  of  C7i{V)  is  contained  in  C7i{V). 

Lemma  3.1  If  two  points  p,r  £V  lie  in  the  same  pocket  AC,-,  then  there  is  a  unique 
external  shortest  path  between  p  and  r  and  it  lies  entirely  in  fCi. 

Proof:  Suppose,  on  the  contrary,  that  some  shortest  path  n  between  p  and  r 
contains  a  point  outside  fCi.  Let  a  (resp.  cj)  be  the  first  (resp.  last)  intersection 
point  of  TT  and  the  lid  q,qi+i.  We  can  obtain  a  shorter  external  path  tt'  by  replacing 
the  portion  of  it  between  c,  and  c/  with  the  segment  cTc/,  contradicting  the  choice 
of  TT.  Hence  tt  =  7r(p,  r)  Hes  entirely  in  AC..  But  it  is  known  that  the  shortest  path 
from  p  to  r  lying  entirely  in  Ki  is  unique  (see  [LP84]).  Thus,  n{p,r)  is  unique. 

D 

The  above  Lemma  shows  that  the  problem  of  finding  an  external  shortest  path 
between  two  points  lying  in  the  same  pocket  is  equivalent  to  finding  the  shortest 
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path  between  two  points  lying  inside  a  simple  polygon.  The  latter  problem  has  been 
extensively  studied  and  can  be  solved  quite  efficiently  [LP84],  [GHLST87]. 

Now  we  focus  our  attention  on  external  shortest  paths  between  points  of  different 
pockets  of  V.  In  the  remainder  of  this  section  we  assume  that  p  e  V  (resp.  r  e  V) 
is  a  point  lying  in  the  pocket  AC,  (resp.  /Cj)  and  i  ^  j . 

Lemma  3.2  The  shortest  left  and  right  paths  from  p  to  r  are  unique.  They  are 
given  by 

TTlip^r)      =      7T{p,qr+i)\\qi^r[q:^\\  ■■■  \\q—[qj\\T^{q:,r)  (3.1) 

■^R{p,r)     =     7r(p,g,)||^75nr||  •••  ||97+2^j+r||7r(gj+i,r),  (3.2) 

where  7r(p,  g,+i)  and  7r(p,  g.)  (resp.  Tr{qj,r)  and  ■!r{qj+i,r))  are  shortest  paths  within 
AC,  (resp.  K.^). 

Proof:  It  is  sufficient  to  show  that  the  path  given  in  the  right-hand  side  of  (3.1) 
is  indeed  strictly  shorter  then  any  other  left  path  from  p  to  r.  As  no  shortest  path 
connecting  p  and  q  meets  the  exterior  of  the  region  enclosed  by  CH{V),  7ri,(p,  r) 
must  at  least  pass  through  all  of  {qi+\,  . . .  ,qj}  or  all  of  {q^+t,  ■  ■  ■  ,Qi}-  However, 
any  external  path  from  p  to  r  that  stays  inside  (or  on)  C'H{V)  but  does  not  pass 
through  some  9,',  i  +  1  <  i'  <  j  is  a  right  path.  Therefore,  7r/,(p,  r)  peisses  through 
each  of  {qi+i,  . . .  ,qj}-  (Being  a  shortest  left  path,  it  visits  each  point  exactly  once.) 
A  similar  argument  shows  that  it  must  visit  qi+i,  ...  ,qj  in  clockwise  order,  hence 
7rL(p, r)  contains  the  portion  of  CT-L{V)  clockwise  between  5,+i  and  ^j,  as  asserted. 
As  for  the  section  of  ■ni{p.,r)  from  p  to  g,+i  (resp.   from  qj  to  r),  it  is  unique  and 

must  stay  completely  in  K-i  (resp.  in  ACj)  by  Lemma  3.1.  The  Lemma  follows. 

D 

To  reiterate,  7rL(p,  r)  (resp.  7rR(p,  r))  consists  of  three  parts: 

(i)  Shortest  path  from  p  to  the  lid  vertex  5,+i  (resp.  g,)  within  AC,-. 

(ii)  The  clockwise  (resp.  counter-clockwise)  portion  of  C'H{V)  between  g,+i  and 
qj  (resp.  q,  and  g^+i). 

(iii)  Shortest  path  inside  AC^  from  the  lid  vertex  g_,  (resp.  qj^\ )  to  r. 

Hence,  7r(p,  r)  also  consists  of  three  parts,  as  above,  because  it  coincides  with  either 
■Ktip.r)  or  ■Knip.r). 

To  compute  the  external  distance  between  two  arbitrary  points  of  V  lying  in 
different  pockets,  we  need  a  fast  procedure  for  determining  the  external  distance 
between  two  points  on  CHiV),  and  a  procedure  for  computing  the  distance  from  a 
lid  vertex  of  a  pocket  AC,  to  an  arbitrary  query  point  of  Ki.  The  external  distance 
between  two  points  on  C'H{V)  can  be  computed  quickly  in  a  straight  forward  man- 
ner by  precomputing  "partial  sums,"  from,  say,  go-    As  for  computing  d{p,q,)  (or 
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d{p,qi+i))  for  p  G  AC,,  in  Lemma  3.1  we  have  shown  that  the  corresponding  a  path 
Hes  entirely  inside  K.,.  Therefore  we  can  use  the  following  result  of  Guibas  et  al. 
[GHLST87]  to  compute  the  distance  from  a  lid  vertex  of  a  pocket  to  a  query  point 
of  that  pocket. 

Given  a  simple  polygon  Q  and  a  vertex  v  of  Q,  the  anchor  of  a  point  x  G  Q(with 
respect  to  v  is  the  last  vertex  on  the  (geodesic)  shortest  path  in  Q  from  i'  to  x,  and 
shortest  path  partition  of  Q  with  respect  to  v  is  the  partition  of  the  boundary  of  Q 
into  maximal  segments  such  that  the  anchor  for  all  points  in  a  segment  remains  the 
same.  Guibas  et  al.  [GHLSTS7]  have  proved 

Theorem  3.3  ([GHLST87])  For  a  given  triangulated  polygon  Q  with  t  vertices 
and  a  vertex  v,  its  shortest  path  partition  with  respect  to  v  can  be  obtained  in  linear 
time.  Once  this  partition  has  been  obtained,  the  geodesic  shortest  distance  from  v 
to  any  point  x  G  Q  can  be  computed  m  time  0(log/).  Moreover,  given  the  anchor 
of  X,  the  geodesic  shortest  distance  can  be  computed  in  0{\)  time. 

Thus,  to  compute  the  external  distance  between  any  two  points  of  V,  first  com- 
pute C'H{V)  to  obtain  the  pockets  of  V\  this  can  be  done  in  0{n)  time  using  the 
algorithm  of  [GYS3],  for  example.  Then  for  each  pocket  AC,,  obtain  the  shortest 
path  partition  of  ACj  from  qi  as  well  as  from  g,+i.  At  this  point,  the  three  portions 
of  di{p,r)  and  of  dR{p,r),  and  hence  d{p,r)  can  be  computed  quickly.  Therefore, 
we  have  shown 

Lemma  3.4  After  linear-time  preprocessing  of  a  triangulated  polygon  V ,  the  exter- 
nal distance  between  any  two  points  of  V  lying  in  different  pockets  can  be  computed 
in  O(logn)  time. 

a 

For  a  point  p  G  AC,,  let  the  left  anchor  (resp.  right  anchor)  denote  the  anchor  of 
p  with  respect  to  q^  (resp.  qs+i).  Observe  that 

(i)  Given  both  left  and  right  anchors  for  x,  y  G  'P,  lying  in  different  pockets, 
d{x,y)  can  be  computed  in  0(1)  time. 

(ii)  The  actual  path  7r(x,  y)  can  be  computed  in  time  0(log  n  +  k),  where  k  is  the 
number  of  turns  in  7r(x,  y). 

Now  we  prove  a  crucial  property  of  external  shortest  paths  which  we  call  mono- 
tonicity  property. 

Lemma  3.5  (Monotonicity  property)  7/r  G  'P{p,s),  then 

dL{p,s)  <  dii{p,s)  =>  dL{p,r)  <  dR{p,r).  (3.3) 
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Figure  3:  Illustrating  the  proof  of  Lemma  3.5:  Four  paths  emanating  from  x 

Proof:  By  Lemma  3.1  the  above  inequality  trivially  holds  if  r  or  5  lie  in  the  same 
pocket  as  p,  therefore  we  only  consider  the  case  when  both  r  and  s  lie  outside  the 
pocket  containing  p.  For  a  contradiction,  suppose  dL{p,s)  <  df{(p,s)  =  di{s,p), 
but  diip^r)  >  dR{p,r)  —  di{r,p).  Thus 

diis.p)  +dL{p,r)  >  diir.p)  +  diip.s).  (3.4) 

By  assumption,  p,  r,  s  occur  in  this  clockwise  order  on  V.  In  particular,  ttl{p-,s) 
intersects  7rx,(r,p).  Let  x  be  point  lying  on  7ri(p,  5)n7r£,(r,p)  but  outside  the  pocket 
of  p;  it  is  easily  seen  that  such  a  point  always  exists.  Denoting  by  tt"^  the  reverse 
of  a  (directed)  path  tt,  7ri,(p,  5)  (resp.  7ri(r,p))  can  be  written  as  7rf^||7r2  (resp. 
7rJ^||7r4)  with  ttj,  ttz,  tts,  1:4  emanating  from  x  (see  Figure  3).  Observe  that  Trf  ^||7r3 
is  a  left  path  from  p  to  r,  hence  ^^(p,  J')  <  iTril  +  Kal.  Similarly,  dL(-s,p)  <  |7r2|  +  |7r4|, 
as  7r^^||7r4  is  a  left  path  from  5  to  p.  Adding  the  two  inequalities,  we  obtain 

dL{p,r)  +  dL{s,p)  <  |7ri|  +  IttjI  +  ^3]  +  |7r4|  =  dL{p,s)  +  (/i,(r,p), 

contradicting  (3.4). 

D 

Observe  that  the  inequalities  in  Lemma  3.5  can  also  be  made  "strict"  on  both 
sides  if  desired. 

Corollary  3.6  For  a  point  p  G  V,  there  exists  a  pair  of  points  p',  p"  G  V  (not 
necessarily  distinct)  such  that 

(t)     \/xeV\p,p'),     dLip,x)     <     dnip^x), 

(it)     Vx6P[p',p"],    dL{p,x)     =     dn{p,x), 

(Hi)     yxeV{p",pl     dL{p,x)     >     dnip.x). 
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Proof:  Letp  e  AC,.  By  Lemma 3.1,  (i/?(p,g,+i)  >  c?l(p,9.+i)  and  (^^(p,^,)  >  dfi{p,qi). 

Since  the  function  A(x)  =  dR{p,x)  —di{p^x)  is  continuous,  the  Intermediate  Value 

Theorem  guarantees  the  existence  of  a  point  r  G  P[g,+i,9,]  such  that  A(r)  =  0,  i.e. 

dL{p,r)  =  dR{p,r).   Let  p'  (resp.  p")  be  the  first  (resp.  last)  point  on  V[qi+-i,qi]  for 

which  A(x)  =  0.    Now,  Corollary  3.6  follows  immediately  from  the  monotonicity 

property. 

D 

Thus,  for  every  point  p  G  "P,  there  exists  a  unique  point  p"  (namely,  the  last 
point  of  V[qi+i,qi]  for  which  A(a:)  =  0)  such  that  for  all  x  G  Vlp^p"],  we  have 
7r(p,  x)  =  nL{p,x),  and  for  all  x  G  'P{p",p],  we  have  7r(p,  x)  =  7r/}(p,  x).  We  call  p" 
the  bifurcation  point  of  p,  and  denote  it  by  /3{p). 

Corollary  3.7  If  y  E  'P\p,l3{p)]  then  for  all  x  G  V[p,y],  dL{x,y)  <  dR{x,y). 

Proof:  Immediate  from  Corollary  3.6  and  the  monotonicity  property. 

D 

Corollary  3.8  Let  p,r  be  two  points  on  an  edge  of  V  such  that  r  lies  after  p  in 
clockwise  direction.  Then,  p,  r,  I3{p),  /9(r)  occur  in  this  clockwise  order  on  V 
(possibly  with  /5(p)  =  /3{r)). 

Proof:  Since  p  and  r  lie  in  the  same  pocket  AC,,  dL{p,r)  <  dfi{p,r)  (and  dL{r,p)  > 
<^r(''iP))-  Thus  I3(p),^(r)  G  'P[r,p].  Ordering  p,  r,  0{r),  /3{p)  contradicts  mono- 
tonicity property  unless  0(p)  =  f3{r).  The  claim  follows. 


4      Computing  external  farthest  neighbors 

In  this  section  we  describe  a  two-phase  algorithm  for  computing  external  farthest 
neighbors  for  all  vertices  of  V.  First  we  determine  a  farthest  neighbor  of  p  among 
the  points  lying  in  the  same  pocket  as  p,  and  then  we  compute  a  member  of  (f)'{p). 
The  farther  of  the  two  will  be  an  external  farthest  neighbor  of  p.  At  the  end  of  this 
section  we  discuss  how  to  modify  the  algorithm  so  that  it  computes  all  external 
farthest  neighbors  for  every  vertex  of  p. 

In  Lemma  3.1  we  have  shown  that  the  external  shortest  path  between  two  points 
lying  in  the  same  pocket  does  not  leave  the  pocket;  this  implies  that  for  every  vertex 
p  G  AC,,  the  problem  of  computing  a  futhest  neighbor  among  the  points  lying  in  the 
same  pocket  as  p  is  equivalent  to  finding  its  farthest  neighbor  inside  the  simple 
polygon  Ki.  Suri  [SuS7],  and  Guibas  and  Hershberger  [GHS7]  have  independently 
given  an  0{tlogt)  algorithm  to  compute  one  internal  geodesic  farthest  neighbor  for 
every  vertex  in  a  simple  polygon  with  t  vertices.  Repeated  application  of  any  of 
these  algorithms  to  each  pocket  of  V  in  turn  yields  the  following  lemma 
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Lemma  4.1   For  every  vertex  p  of  V,   a  farthest  neighbor  of  p  among  the  points 
lying  m  its  pocket  can  be  computed  in  time  0(n  log  n). 

D 

In  the  remainder  of  this  section  we  develop  an  efficient  algorithm  to  compute  a 
representative  of  (^'(p)  for  each  vertex  p  ^  V.  For  a  point  p  G  X^,,  let  l(p)  =  g,+i 
and  r{p)  =  q,.  Let  C{p)  (resp.  7^(p))  denote  the  vertices  of  V  lying  in  V[l{p), /3{p)] 
(resp.  V[0{p),r(p)])  ordered  in  clockwise  direction.  Let  6l(p)  denote  the  last  vertex 
p'  G  C{p)  such  that 

dL{p,p')  =  max  {dL{p,x)\  X  e  C{p)}  . 

Similarly,  let  <5/}(p)  denote  the  first  vertex  p"  G  "^(p)  such  that 

dniP^p")  =  max  {dR{p,x)\x  G  "^(p))  . 

Intuitively,  6l{p)  (resp.  <5/?(p))  represents  a  vertex  of  V  farthest  from  p  among 
the  vertices  lying  in  'P[p,l3{p)]  (resp.  V[j3{p),p\). 

Lemma  4.2   For  any  point  p  E  V,  at  least  one  of  8l{p),  l^ip),  8r{p)  is  in  (f>'{p). 

Proof:  Suppose  x  G  (f)'{p)nV[l{p),  f3{p)].  By  definitions  of  <j)'{p)  and  ^(p),  we  must 
have 

dL{p,x)  =  sup{dL(p,y)|y  G  V[l{p), 0{p)]]. 

But  it  follows  from  Lemma  1  of  [PSRS6]  that  d{p,y)  —  di{p,y)  is  a  convex  func- 
tion of  y  on  any  straight-line  segment  of  'P[/(p),/3(p)],  and  hence  d{p,y)  attains  its 
maximum  at  either  of  the  endpoints.  Thus  x  G  C{p)  U  {/?(p)},  as  asserted.  The  case 
X  G  <f>'{p)  n  'P[(3{p),r{p)]  is  similar. 

D 

Observe  that  the  proof  of  the  previous  Lemma  implies  that  the  point(s)  of  {<5l(p), 
/9(p),  <5r(p)}  farthest  from  p  must  lie  in  </i>'(p).  Therefore,  a  representative  of  <f>'{p) 
can  be  located  by  computing  Si{p),  I3{p),  and  Sji{p). 

Assume  that  V  is  given,  CTiiV)  has  been  computed,  and  the  preprocessing  for 
external  distance  queries  has  been  performed.  We  will  first  describe  how  to  compute 
/3(p)  for  all  vertices  p  G  'P  and  then  present  an  algorithm  that  determines  6l{p)  for 
each  vertex  p.  The  procedure  for  computing  6ji{p)  is  similar. 

4.1      Locating  bifurcation  points 

Let  n|  (resp.  H^)  denote  the  set  of  endpoints  of  the  segments  of  the  shortest-path 
partition  of  AC,  from  g,  (resp.  q,+i).  Let  Ui  —  po,  and  U2,  ...  ^Un  be  the  points  of 
n  =  [J(n|  U  n[),  sorted  along  V  in  the  clockwise  direction.  (Note  that  the  set  V  of 
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vertices  of  V  is  contained  in  IT  by  construction.)  11  partitions  V  into  0{n)  segments 
such  that  the  left  and  right  anchors  of  a  point  x  remain  the  same  as  it  varies  over 
a  single  segment  (cf.  [GHLST87]).  If  v  is  the  right  anchor  of  x,  then  dii{p,x)  = 
dR{p,v)  +  d{v,x),  where  d{p,v)  is  fixed  as  x  ranges  over  utUjt+i,  and  {d{v,x)y 
depends  quadratically  on  the  position  of  x  on  tZpuT+T-  The  function  di{x,y)  has  a 
similar  form. 

Lemma  4.3  It  is  possible  to  compute  H  and  determine  left  and  right  anchors  for 
each  point  of  H  and  for  each  segm,ent  of  V  —  H  in  triangulation  time. 

Proof:  n|  and  !![  can  be  computed  in  desired  time  as  in  Lemma  3.4.  IT  is  deter- 
mined by  merging  fl'  and  H^  for  each  i  and  concatenating  the  resulting  sequences. 
As  left  (resp.  right)  anchor  of  a  point  p  €  HJ"  (resp.  n|)  depends  only  on  the  relative 
position  of  p  among  points  of  n|  (resp.  II^),  for  each  point  on  11  both  anchors  can 

be  easily  determined  at  merge  time.  Segments  oi  V  —  U.  can  be  handled  similarly. 

D 

For  X,  y  ^  V,  let  A(x,y)  =  d/{(a-,y)  —  dL{x,y).  By  Corollary  3.6  the  sign  of 
A(p, z)  changes  monotonically  as  z  moves  along  V[l{p),r{p)],  therefore  there  is  a 
unique  u,  G  11  such  that  A(p,  u,)  >  0  and  A(p,  Ui+i)  <  0.  Since  0{p)  is  the  last  point 
z  G  'P[/(p),^(p)]  such  that  A{p,z)  =  0,  /?(/>)  G  u,u,+i.  Moreover,  given  i,  0{p)  can 
be  determined  in  constant  time,  as  the  precise  analytic  form  of  A(p,  z)  on  u,u,+i  is 
known  (refer  to  the  discussion  above).  Thus  we  can  compute  bifurcation  points  for 
all  vertices  of  T'  as  follows:  We  first  compute  /5(po)  by  scanning  the  list  11  starting 
from  u,o  =  /(po)  until  the  first  i  such  that  A(po,w,)  <  0  and  then  locating  ^(po)  in 
the  segment  u,_iu,.  Now  inductively  assume  that  we  have  computed  /3(pj_i)  and  it 
lies  in  uiul+i-  By  Corollary  3.8,  ^(pj)  G  T'[/?(pj_i),Pj_i],  therefore  we  scan  11  from 
u^,  find  the  first  i  such  that  A(pj,u,)  <  0,  and  locate  /?(pj)  in  ui~[Ui. 

Lemma  4.4  The  above  algorithm  computes  0{p,)  for  all  vertices  pi  E  V  in  linear 
time. 

Proof:  Correctness  of  the  algorithm  follows  immediately  from  Corollary  3.8  and 
from  the  previous  discussion.  Observe  that  each  external  distance  query  used  in  the 
algorithm  is  between  two  points  in  different  pockets  and  requires  0(1)  time;  this  is 
because  the  points  lie  in  IT  and  we  have  assumed  that  left  and  right  anchors  for  all 
such  points  have  been  precomputed.  Therefore  the  running  time  of  the  algorithm 
is  proportional  to  the  number  of  times  A  is  evaluated.  Corollary  3.8  implies  that 
every  evaluation  of  A  is  associated  with  advancing  either  in  the  list  of  vertices  or 
in  the  fist  IT  of  candidate  points.  Since  we  traverse  the  list  IT  at  most  twice  (once 
to  compute  ^(po)  and  at  most  once  to  compute  ^(pi),  ...  ,/3(p„_i),  as  /3(p„_i)  G 
'P{Po,0{po)]),  A  is  evaluated  at  most  0{n)  times.  Hence  the  total  time  spent  in 
computing  ^(p.)  for  all  vertices  of  V  is  at  most  0{n). 

D 
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Figure  4:   Left  distance  from  x  to  vertices  pt  G  V^i^jiPj];  bullets  denote  5[9,,Pj]  C 

4.2      Computing  ^l  and  8^ 

In  this  section  we  describe  how  to  compute  8l  for  all  vertices  of  V\  Sr  is  computed 
similarly.  For  p,  g  G  "P,  let  V[p,  9]  denote  the  sequence  of  vertices  lying  on  P\p,q\  in 
clockwise  direction. 

Recall  that  Sl{p)  is  the  last  vertex  of  C.{p)  farthest  from  p.  A  brute-force  al- 
gorithm for  computing  ^^(p,)  for  a  vertex  pi  of  V  would  determine  the  external 
distance  from  p,  to  every  point  in  £(pi),  and  would  choose  the  one  farthest  from 
p,.  However,  the  worst  case  time  complexity  of  computing  8i{pi)  for  all  vertices  of 
V  using  this  approach  is  quadratic.  Consequently,  we  now  describe  the  intuition 
behind  the  changes  that  will  transform  this  straight  forward  procedure  into  an  ef- 
ficient algorithm.  We  regard  C{pi)  as  a  queue  of  points  ordered  along  V,  with  its 
front  at  /(p,)  and  rear  towards  I3{pi).  Observe  that,  as  one  advances  from  pi  to 
Pi+i,  £(p,+i)  can  be  obtained  from  £(p,)  by  deleting  V[l{pi),l{pi^i))  from  its  front 
and  appending  F(^(p,), /?(p,+i)]  to  the  rear.  (The  former  step  is  non- vacuous  only 
if  pi+i  is  a  lid  vertex  because  otherwise  /(p,)  =  /(pi+i).)  However,  this  modifica- 
tion alone  does  not  improve  the  performance  of  the  algorithm,  since  for  each  p,  it 
is  still  necessary  to  examine  every  point  in  C(pi)  when  looking  for  a  maximum  of 
the  function  di{pi,x).  The  next  improvement  follows  from  the  observation  that 
only  some  specific  points  of  £(p,)  have  a  chance  of  being  8L{pi).  More  precisely, 
we  define  a  subsequence  of  £(p,)  so  that  8i{pi)  is  guaranteed  to  be  the  first  vertex 
of  this  subsequence.  In  addition,  the  corresponding  subsequence  of  £(p,+i)  can  be 
built  from  that  of  C{pi)  by  deleting  from  its  front  all  vertices  (if  any)  that  lie  in 
V[/(p,), /(p,+j)),  and  appending  to  the  rear  selected  vertices  of  y(^(p,),/?(pi+i)]. 
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Figure  5:  Illustration  of  Lemma  4.5 

More  formally,  let  ^^[9i,Pj]  be  the  sequence  of  consecutive  vertices  of  7^,  starting 
at  a  vertex  g,  lying  on  CT-l{V),  and  ending  at  pj.  Then  the  peak  sequence  of  V[qi,pj], 
denoted  by  S[5,, Pj],  consists  of  cill  vertices  pk  G  ^^[^i, Pj]  such  that 


Vpt  G  V{pk,pj],  diix.pk)  >  dL{x,pt), 


(4.1) 


where  x  is  an  arbitrary  point  of  V[pj,qi]  (see  Figure  4).  The  peak  sequence  S[qi,pj] 
is  well-defined  in  the  sense  that  the  condition  (4.1)  is  independent  of  the  choice  of 
X  G  'P\pj,Qi],  as  the  following  lemma  shows  (taking  Zi  =  p^,  22  =  p,  and  y  to  be  an 
alternative  choice  of  x): 

Lemma  4.5  Let  zi  and  22  be  two  points  on  V,  and  let  x,  y  £  V[zi,  22]  be  two  points 
such  that  V[x,  22]  Pi  V[y,  22]  contains  a  vertex  qi  €  CH{V)  (see  Figure  5).    Then 

dLix,Z2)  >  dL{x,zi)  -^  di{y,Z2)  >  diiy.Zi). 

Proof:  Since  9,  e  V[x,  22]  n  Cn{V),  q,  lies  on  7r^(a;,  22),  and  so  diix,  22)  =  ^^(x,  g.) 
+  dLiq^,Z2).  Similarly,  diiy^z^)  =  di{y,qi)  +  diiq^.z^),  dL{x,Zi)  =  diix.qi)  + 
dL{qi-,zi),  and  dL{y,zi)  =  dL{y,qi)  +  dL{q„Zi).  Therefore 

dL{x,Z2)>  diix.zi)      <»     dLix,q,)  +  dL{q,,Z2)>  dL{x,q,)  +  dLiqi,Zi) 

^     dLiy,q^)  +  dL{q„Z2)  >  dciy^q,)  +  dL(q„Zi) 
^     dL{y,Z2)>  diiy^Zi). 

a 

Fixing  a  point  x  €  V[p.j,q,],  observe  that  di{x,p)  strictly  decreases  as  p  moves 
clockwise  through  the  elements  of  S[g,,pj].  Thus  we  have 
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Input:       V  preprocessed  for  externa]  distance  queries, 
and  /3(pi)  for  all  p,  G  V . 
Output:  Si{p,)  for  each  pi  G  V. 

Q  ■■=  (Po);    X  :=  Po 
for  ji  =  0  . . .  n  —  1  do 

forpk  eV{x,0{Pj)]  do 

while  Q^<i)k.  dL{Pj,rear{Q))  <  di{pj,pk)  do 

delete. Tear{Q) 
end  while 
insertjrear{  Q,  p^.) 

end  for 

while  Q^4>k  frontiQ)  G  V[/(pj_i),/(pj))  do 

delete.front{Q) 
end  while 

^l(Pj)  :=  front{Q)-     x  =  fi{Pj) 
end  for 


Figure  6:  Computing  StiPi)  for  vertices  of  V 

Lemma  4.6    The  first  element  of  the  peak  sequence  of  C{p,)  is  8i{pi). 

D 

Hence  <5l(p,)  can  be  determined  by  computing  5[/(p,),/?(p,)],  so  we  only  have  to 
show  how  to  calculate  S[/(p,),;3(p,)],  for  every  vertex  p,  of  V.  Given  a  vertex  po,  a 
clockwise  scan  of  V [/(po),  /?(Po)]  is  sufficient  to  compute  5[po,  /9(po)]  and  thus  (5l(po). 
The  following  two  lemmas  suggest  an  efficient  method  for  computing  S[/(p,), /9(p,)] 
from  5[/(p._i),/3(pi_i)]. 

Lemma  4.7  If  V[qt,pj]  C  Vlq^-^pj],  then  S[q,i^pj\  is  simply  Slq^^pj]  truncated  to 
within  V[qii,pj\. 


Proof:  Immediate  from  the  definition  of  a  peak  sequence. 


D 


Lemma  4.8  The  peak  sequence  Sj  =  5[5,,pj+i]  can  he  obtained  from  S\  =  S[qi,pj] 
by  first  deleting  all  points  pk  of  Si  for  which  d[^{qi,pk)  <  di(qi,pj^i),  and  then 
inserting  Pj+i- 
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Proof:  Observe  that  pj+i  6  52  and  S2  —  {Pj+i}  Q  Si.  Suppose  pk  6  ^[^i, Pj]  is  an 

element  of  5].  Then  diiqi^pk)  >  dL{qi,pk')  whenever  k  <  k'  <  j .  Thus  pk  is  also  52 

if  and  only  if  diiqiyPk)  >  dL{qiTPj+i)i  as  asserted.  The  choice  of  5,  as  the  reference 

point  in  'P[pj+i,g,]  was  purely  arbitrary. 

D 

Note  that  the  elements  deleted  from  Si  to  obtain  52  form  a  suffix  of  Si . 

To  compute  5[/(pj),^(pj)]  from  5[/(p,_i), /?(p,_i)],  first  extend  this  peak  se- 
quence to  S[l{pj-i),  j3{p.j)]  by  repeatedly  applying  Lemma  4.8  to  the  vertices  of 
V{^{pj-i),  (3{p.j)],  and  then  truncate  5[/(pj_i),/3(pj)]  to  S[l{pj),  I3{pj)]  as  in  Lemma 4.7. 
A  more  formal  description  of  this  algorithm  is  given  in  Figure  6. 

The  correctness  of  the  algorithm  follows  from  the  above  discussion.  As  for  the 
running  time  of  the  algorithm,  we  have: 

Lemma  4.9  The  time  required  to  compute  Siipi)  for  all  vertices  pi  G  "P,  after 
preprocessing  V  for  external  distance  queries,  is  at  most  0{n). 

Proof:  It  follows  from  Lemma  4.4  that  the  bifurcation  points  of  all  vertices  of  V  can 
be  calculated  in  0{n)  time.  As  for  updating  the  queue  Q,  V  is  scanned  at  most  twice 
(once  for  computing  Si{pq),  and  once  for  computing  8l  for  the  remaining  vertices), 
which  implies  that  onlj^  0(n)  points  are  inserted  in  Q.  Since  each  element  is  deleted 
from  the  queue  only  once,  only  0{n)  queue  insertions  and  deletions  are  performed 
and  each  requires  0(1)  time.  Moreover,  each  distance  query  is  immediately  followed 
by  a  queue  insertion  or  deletion  or  advance  in  the  list  of  vertices,  hence  all  Sl  can 
be  computed  in  0{n)  time  after  the  initial  preprocessing. 

D 

Using  the  same  method  we  can  compute  6r  for  all  vertices  in  0(n)  time.  Therefore, 
we  conclude  that 

Theorem  4.10  One  point  of  4>'{Px)  for  each  vertex  p,  G  V  can  be  computed  in 
triangulation  time. 

a 

Remark:  The  above  algorithm  computes  only  a  representative  of  <f>'{p,).  However, 
it  can  be  easily  extended  so  that  it  returns  all  elements  of  (l)'{p,)  in  time  0{n  +  k'), 
where  k'  =  J2p,ev  \4>'{p,)\-  First  we  run  the  above  algorithm  to  determine  which  of 
{^l(p.),  0{Pi),  8r(p,)}  is  farthest  from  p,.  If  6c{pi)  (resp.  6R{pi))  is  farthest,  then  we 
compute  all  vertices  in  V[l(p,),  0{p,)]  (resp.  V[0{p,),r{p,)])  which  are  farthest  from 
p,  by  running  a  slightly  modified  version  of  our  algorithm  that  uses  ">"  instead 
of  ">"  in  the  defining  property  (4.1)  of  peak  sequences.  In  this  variation  of  the 
algorithm,  the  set  of  all  vertices  of  £(p.)  farthest  from  p,  forms  a  contiguous  initial 
portion  of  Q  and  thus  can  be  extracted  from  it  (if  necessary)  in  time  proportional 
to  the  number  of  such  points. 
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Figure  7:  A  simple  polygon  V  with  k'  =  Q.{n^) 

It  can  be  easily  shown  that  in  the  worst  case  k'  =  fi(n^).  For  the  lower  bound, 
consider  the  simple  polygon  V  shown  in  Figure  7.  The  points  Oi,  ...  ,  Ot,  for  t  —  |— 5 
lie  on  a  circular  arc  with  center  c/  and  radius  p,  and  the  points  6i,  ...  ,  6(  lie  on 
another  circular  arc  with  center  Cy  and  radius  p.  If  we  choose  di  and  dr  so  that 
cidi  =  Cydr.,  diici^dr)  =  dfi{di^Cr),  and  di  (resp.  dr)  is  not  visible  from  at  (resp.  6i), 
then  it  is  easy  to  see  that  for  any  i  <  t,  4>{a,)  —  {bi,  ...  ,fe()  and  (resp.  (f>(bt)  = 
{ai,...  ,a(}.  Hence  X^p.e-P  I  </"'(?.)  I  =  f^(n^). 

Corollary  4.11    The  total  time  required  to  compute  one  external  farthest  neighbor 
4>{pi)  for  each  vertex  pi  G  V  is  at  most  O{n\ogn). 


Remark:  Note  that  the  only  part  of  the  algorithm  that  required  0{n  log  n)  time 
is  the  computation  of  a  farthest  neighbor  of  p,  among  the  points  lying  in  the  same 
pocket  as  p,,  for  each  p,  E  V,  as  the  preprocessing  of  our  algorithm  requires  only 
triangulation  plus  linear  time  and  after  preprocessing  a  representative  of  4>'{p,)  for 
each  Px  EV  can  be  computed  in  0{n)  time.  Therefore,  a  more  efficient  algorithm  to 
compute  the  geodesic  farthest  neighbors  would  immediately  yield  a  faster  algorithm 
for  the  external  farthest  neighbors.  Moreover,  we  can  modify  the  algorithm  as 
discussed  above  together  with  a  variant  of  the  algorithm  by  Suri  [Su87]  or  that  given 
by  Guibas  and  Hershberger  [GH87]  to  yield  an  0{k  +  n  log  n)  algorithm  for  obtaining 
all  external  farthest  neighbors  for  each  vertex  p  G  V,  where  k  =  X!!pev  l?^!/^)!- 
As  for  computing  the  external  diameter  of  7^,  it  follows  from  Lemma  2.1  that 

Theorem  4.12    The  external  diamater  of  a  simple  n-gon  can  be  computed  in  time 
0(n  log  n)  using  0{n)  space. 
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5      Conclusions 

In  this  paper,  we  provided  an  0(n  log  n)  algorithm  for  finding  an  external  geodesic 
farthest  neighbor  for  every  vertex  of  a  given  simple  n-gon.  We  also  showed  that 
this  algorithm  can  be  modified  to  output  all  external  geodesic  farthest-neighbors 
for  all  the  vertices  of  the  polygon,  and  that  it  can  be  used  to  compute  the  external 
diameter  of  the  polygon  in  0(n  log  n)  time.  For  finding  the  external  diameter, 
the  time  complexity  of  our  algorithm  improves  upon  that  given  by  Samuel  and 
Toussaint  [ST87]  by  an  0{n/\ogn)  factor.  However,  the  best  known  lower  bound 
for  computing  the  external  geodesic  diameter  is  the  trivial  0{n)^  and  obtaining  an 
optimal  bound  remains  tantalizingly  open. 

A  careful  analysis  of  our  algorithm  shows  that  if  there  is  an  algorithm  that 
computes  the  internal  geodesic  fcirthest  neighbors  of  the  vertices  of  a  simple  n-gon 
in  T{n)  time,  then  our  algorithm  can  use  it  as  a  subroutine  to  compute  one  external 
farthest  neighbor  of  every  vertex  of  a  triangulated  simple  n-gon  in  time  0{T{n)  +  n) 
—  0{T(n)).  Consequently,  an  improvement  in  the  time  complexity  of  finding  the 
internal  geodesic  farthest  neighbors  will  yield  the  corresponding  improvement  in  the 
running  time  of  finding  the  external  farthest  neighbors. 
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